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Abstract 

Recently the possibility of generating nonclassical polariton states by means of parametric scat- 
tering has been demonstrated. Excitonic polaritons propagate in a complex interacting environment 
and contain real electronic excitations subject to scattering events and noise affecting quantum 
coherence and entanglement. Here we present a general theoretical framework for the realistic in- 
vestigation of polariton quantum correlations in the presence of coherent and incoherent interaction 
processes. The proposed theoretical approach is based on the nonequilibrium quantum Langevin 
approach for open systems applied to interacting-electron complexes described within the dynamics 
controlled truncation scheme. It provides an easy recipe to calculate multi-time correlation func- 
tions which are key-quantities in quantum optics. As a first application, we analyze the build-up 
of polariton parametric emission in semiconductor microcavities including the influence of noise 
originating from phonon induced scattering. 
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I. INTRODUCTION 



Entanglement is one of the key features of quantum information and communication 
technology Parametric down-conversion is the most frequently used method to gener- 
ate highly entangled pairs of photons for quantum-optics applications, such as quantum 
cryptography and quantum teleportation. Rapid development in the field of quantum infor- 
mation requires monolithic, compact sources of nonclassical photon states enabling efficient 
coupling into optical fibres and possibly electrical injection. Semiconductor-based sources 
of entangled photons would therefore be advantageous for practical quantum technologies. 
Moreover semiconductors can be structured on a nanometer scale, and thus one may produce 
materials with tailored properties realizing a wide variety of physically distinct situations. 
However semiconductor heterostructures constitute a complex interacting environment in- 
volving charge, spin, and lattice degrees of freedom, hence suited to serve as prototype 
systems where quantum-mechanical properties of many interacting particles far away from 
equilibrium can be studied in a controlled fashiony]. It has been demonstrated that very 
large x^^^ resonant polaritonic nonlinearities in wide-gap semiconductors and in semiconduc- 
tor microcavities can be used to achieve parametric emission 3|, |J] . 

Polaritons are mixed quasiparticles resulting from the strongly coupled propagation of 
light and collective electronic excitations (excitons) in semiconductor crystals. Although 
spontaneous parametric processes involving polaritons in bulk semiconductors have been 
known for decades j^, the possibility of generating entangled photons by these processes was 
theoretically pointed out only lately [sl. This result was based on a microscopic quantum 
theory of the nonlinear optical response of interacting electron systems relying on Jthe dy- 
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namics controlled truncation scheme p] extended to include light quantization 



m. The 



above theoretical framework was also applied to the analysis of polariton parametric emis- 
sion in semiconductor microcavities (SMCs)!?, 9]. A SMC is a photonic structure designed 
to enhance light-matter interactions. The strong light-matter interaction in these systems 
gives rise to cavity polaritons which are hybrid quasiparticles consisting of a superposition of 
cavity photons and quantum well excitons [lO]. Demonstrations of parametric amplification 
and parametric emission in SMCsP 



, together with the possibility of ultrafast optical 



manipulation and ease of integration of these microdevices, have increased the interest on 

, llSl . Il6| . In 2004, ex- 
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perimental evidence for the generation of ultraviolet polarization-entangled photon pairs by 
means of biexciton resonant parametric emission in a single crystal of semiconductor CuCl 
has been reported 17|. Short-wavelength entangled photons are desirable for a number of 
applications as generation of further entanglement between three or four photons. In 2005 
an experiment probing quantum correlations of (parametrically emitted) cavity polaritons 
by exploiting quantum complementarity has been proposed and realized[l^. Specifically, 
it has been shown that polaritons in two distinct idler modes interfere if and only if they 
share the same signal mode so that which-way information cannot be gathered, according to 
Bohr's quantum complementarity principle. In 2006 a promising low-threshold parametric 
oscillation in vertical triple SMCs with signal, pump and idler waves propagating along the 



vertical direction of the nanostructure has been demonstrated 

The crucial role of many-particle Coulomb correlations in semiconductors marks a pro- 
found difference from dilute atomic systems, where the optical response is well described by 
independent transitions between atomic levels, and the nonlinear dynamics is governed only 
by saturation effects due to the Pauli exclusion principle. In planar SMCs, thanks to their 
mutual Coulomb interaction, pump polaritons generated by resonant optical pumping may 
scatter into pairs of polaritons (signal and idler) 0, [s], [3] , they are determined by the two cus- 
tomary energy and wave vector conservation conditions 2kp = k^ -|- kj and 2£'kp = -f'ks + -^ki 
depicting an eight-shaped curve in momentum space. At low pump intensities they are ex- 
pected to undergo a spontaneous parametnc ptoeess dnven by vaenunj-flnetnation, whereas 
at moderate intensities they display self-stimulation and oscillation [4]. However they are 
real electronic excitations propagating in a complex interacting environment. Owing to the 
relevance of polariton interactions, and also owing to their interest for exploring quantum 
optical phenomena in such a complex environment, theoretical approaches able to model ac- 
curately polariton dynamics including light quantization, losses and environment interactions 
are highly desired. The analysis of nonclassical correlations in semiconductors constitutes a 
challenging problem, where the physics of interacting electrons must be added to quantum 
optics and should include properly the effects of energy relaxation, dephasing, and noise, 
induced by electron-phonon interaction [2^. 

Previous descriptions of polariton parametric processes make deeply use of the picture 
of polaritons as interacting bosons. These theories have been used to investigate paramet- 
ric amplifications, parametric luminescence, coherent control, entanglement and parametric 
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scattering in momentum space {12!. Il4l . Il5l . Il9l . |2l|. 



It is worth noting that in a reahstic environment phase- coherent nonhnear optical pro- 
cesses involving real excitations compete with incoherent scattering as evidenced by exper- 
imental results. In experiments dealing with parametric emission, what really dominates 
emission at low pump intensities is the photoluminescence (PL) due to the incoherent dy- 
namics of the single scattering events driven by the pump itself and the Rayleigh scattering 
of the pump due to the unavoidable presence of structural disorder. The latter process is 
elastic and can thus be spectrally filtered in principle, moreover it is confined in k-space to 
a ring of in-plane wave vectors with almost the same modulus of the pump wave vector. On 
the contrary PL, being not an elastic process, cannot be easily separated from parametric 
emission. Only once the pumping become sufficient the parametric processes start to reveal 
themselves and to take over pump-induced PL as well. Indeed, usually, parametric emission 
and standard purnp-induced PL cohabit as shown by experiments at low and intermediate 
excitation density 12] • Moreover, in order to address quantum coherence properties and 
entanglement 171] the preferred experimental situations are those of few-particle regimes, 
namely coincidence detection in photon counting. In this regime, the presence of incoherent 
noise due to pump-induced PL tends to spoil the system of its coherence properties lowering 
the degree of nonclassical correlations. The detrimental influence of incoherent effects on the 
quantum coherence properties is also well evidenced in the measured time-resolved visibility 
shown in Ref. Il6|. At initial times visibility is suppressed until parametric emission prevails. 
Thus, a microscopic analysis able to account for parametric emission and pump-induced 
PL on an equal footing is highly desirable in order to make quantitative comparison with 
measurements and propose future experiments. Furthermore a quantitative theory would 
be of paramount importance for a deeper understanding of quantum correlations in such 
structures aiming at seeking and limiting all unwanted detrimental contributions. 

The dynamics controlled truncation scheme (DOTS) provides a (widely adopted) start- 
ingj)oint for the microscopic theory of light-matter interaction effects beyond mean-field 







2|, |6j , supplying a consistent and precise way to stop the infinite hierarchy of higher-order 
correlations which always appears in the microscopic approaches of manybody interacting 
systems. In 1996 the DOTS was extended in order to include in the description the quan- 
tization of the electromagnetic field This extension has been applied to the study of 



quantum optical phenomena in semiconductors as polariton entanglement 



However, in 
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these works damping has been considered only at a phenomenological leveL 

In this paper we shall present a novel approach based on a DCTS-nonequilibrium quan- 
tum Langevin description of the open system in interaction with its surroundings. This 
approach enables us to include on an equal footing the microscopic description of the scat- 
tering channels competing with the coherent parametric phenomena the optical pump in- 
duces. We shall apply our method in order to perform a more realistic description of light 
emission taking into account nonlinear parametric interactions, light quantization, cavity 
losses and polariton-phonon interaction. The developed theoretical framework can be natu- 
rally extended to include other incoherent scattering mechanisms such as the interaction of 
polaritons with thermal free electrons 37|] . As a first application of the proposed theoretical 
scheme, we have analyzed the time-resolved and time-integrated build-up of polariton para- 
metric emission in semiconductor microcavities including the influence of noise originating 
from phonon induced scattering. The presented numerical results clearly evidence the role 
of incoherent scattering in parametric photoluminescence and thus show the importance of a 
proper microscopic analysis able to account for parametric emission and pump-induced PL 
on an equal footing. We also exploit the present approach to calculate the emission spectra 
as a function of the pump power density. The spectra display a significant line-narrowing 
as well as parametric emission starts to prevail. 

The paper is organized as follows. In Sec [Tll starting from a DCTS theory for semi- 
conductor microcavities 22|], we present a theory of x^^^ optical nonlinearities in terms of 
interacting polaritons. The latter focuses mainly on the nonlinear part in order to model 
coherent optical parametric processes and the damping is included only phenomenologically. 
In Section IIIII we apply a nonequilibrium quantum Langevin treatment of damping and 
fluctuations in an open system, originally proposed by Lax. Section [IV] will be devoted to 
the microscopic calculation of phonon-induced scattering rates and polariton PL within a 
second order Born-Markov approximation. In Sec. |V]we shall present a quantum Langevin 
description of parametric emission including incoherent effects; particular attention will be 
devoted to the case of single pump feed, whose results will be the subject of Sec. IVD Finally 
in Sec. IVIII we shall summarize and draw some conclusions. 
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II. DYNAMICS CONTROLLED TRUNCATION SCHEME FOR INTERACTING 
POLARITONS 



The system under investigation consists of one or more (uncoupled) QWs grown inside a 
semiconductor planar Fabry-Perot resonator. For the quasi-2D interacting electron system 
we adopt the usual semiconductor model Hamiltonian [2] which can be expressed as 

i^e = ^ ^Na I Na){Na I , (1) 

Na 

where the eigenstates of He have been labeled according to the number N of electron-hole 
{eh) pairs. The state | = 0) is the electronic ground state, the = 1 subspace is the 
exciton subspace with the additional collective quantum number a denoting the exciton 
energy level n and the in-plane wave vector k. The set of states with N = 2 determines the 
biexciton subspace. We treat the planar-cavity field within the quasimode approximation, 
the cavity field is quantized as though the mirrors were perfect: 

Hc = '^ hukaldk , (2) 

k 

and the resulting discrete modes are then coupled to the external continuum of modes by 
an effective Hamiltonian 

Hp = ^hte, - 4k «k) , (3) 

i=i,2,k 

where j labels the two mirrors and tj determines the fraction of the field amplitude passing 
the cavity mirror, E^^^ (4^"*) positive (negative) frequency part of the coherent input 

light field. The coupling of the electron system to the cavity modes is given within the usual 
rotating wave approximation 

Hi = - Y,^Vr,^,^iB^^^ + H.c. (4) 

nk 

5„k is the exciton destruction operator and can be expanded as well in terms of the energy 
eigenstates of the electron system. For later convenience, the exciton and photon operators 
are normalized so that B\.B]^ and aj^fik are operators corresponding to the number of particles 
within a Bohr-radius two-dimensional disk (vra^) at a given k. 

We start from the Heisenberg equations of motion for the exciton and photon operators. 
In the DOTS spirit, we keep only those terms providing the lowest nonlinear response (x''^^) 
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in the input light field [22| . We assume the pump polaritons driven by a quite strong coherent 
input field E^^ = {e[^J) consisting of a classical (C-number) field, resonantly exciting the 
structure at a given energy and wave vector, kp. We are interested in studying polaritonic 
effects in SMCs where the optical response involves mainly excitons belonging to the IS 
band with wave vectors close to normal incidence, Ikl ^ — . We retain only those terms 

' ' ' fix 

containing the semiclassical pump amplitude twice, thus focusing on the "direct" pump- 
induced nonlinear parametric interaction. One ends up with a set of coupled equations 
of motion exact to the third order in the exciting field. While a systematic treatment of 
higher-order optical nonlinearities would require an extension of the equations of motion, 
a restricted class of higher-order effects can be obtained from solving these equations self- 
consistently up to arbitrary order as it is usually employed in standard nonlinear optics. 
This can be simply accomplished by replacing, in the nonlinear sources, the linear excitonic 
polarization and and light field operators with the total field. From now on, that the 
pump-driven terms (e.g. the B and a at kp) are C-numbers coherent amplitudes like the 
semiclassical electromagnetic pump field, we will make such distinction in marking with a 



'hat" the operators only. It yields [22 1 



5k = -iuj^Bk-isk + iVak-iR^^ , (5a) 
= -iuj^cik + iVl3k + tc^k" ; (5b) 

where {i = x, c) are the energies of QWs excitons and cavity photons. The intracavity 
and the exciton field of a given mode k are coupled by the exciton-cavity photon coupling 
rate V. The relevant non-linear source term, able to couple waves with different in-plane 
wave vector k, is given by R^^ = {R^* + R^)/Neff; where the first term originates from 
the phase-space filling of the exciton transition, 

= ^B^^a^^Bl ; (6) 

''sat 

being Usat = 7/16 the exciton saturation density and kj = 2kp — k. Ne/f depends on the 
number of wells inside the cavity and their spatial overlap with the cavity-mode. Inserting a 
large number of QWs into the cavity results also in increasing the photon-exciton coupling 



rate V = V\ ^jN^ff^ where Vi is the exciton-photon coupling for 1 QW. R^ is the Coulomb 
interaction term. It dominates the coherent xx coupling and for co-circularly polarized waves 
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(the only case here addressed) can be written as 

-i f dt'F{t-t')B^^{t')B^^(t')\ (7) 

where V^x — QEb/Tc, being Eh the exciton binding energy. Equation ([7]) includes the in- 
stantaneous mean-field xx interaction term and a non-instantaneous term originating from 
four-particle correlations. These equations show a close analogy to those derived in [3], 
addressing the bulk case. In addition to that former result, in the present formulation we 
succeed in dividing rigorously (in the DCTS spirit) the Coulomb-induced correlations into 
mean-field and four-particle correlation terms. Moreover the pump-induced shift due to 
parametric scattering reads 

-2iBl^{t) f dt'F{t-t')B^(t')B^^{t'). (8) 

J — oo 

Equation <^ can be written in compact form as 

i3k = -^^^;^^i3k + 4"-^7^i^^; (9) 

( B^.\ ( -v\ . ( \ ( R^A 

where Bk= \ \ , fll^ = \ , ^k" = h and 7^^^ = \ ^ \. When 

the coupling rate V exceeds the decay rate of the exciton coherence and of the cavity field, 
the system enters the strong coupling regime. In this regime, the continuous exchange of 
energy before decay significantly alters the dynamics and hence the resulting resonances of 
the coupled system with respect to those of bare excitons and cavity photons. As a conse- 
quence, cavity-polaritons arise as the two-dimensional eigenstates of The coupling rate 
V determines the splitting (~ 2V) between the two polariton energy bands. This nonper- 
turbative dynamics including the interactions (induced by R^^) between different polariton 
modes can be accurately described by Eq. ([5]). Nevertheless there can be reasons to prefer a 
change of bases from excitons and photons to the eigenstates of the coupled system, namely 
polaritons. An interesting one is that the resulting equations may provide a more intu- 
itive description of nonlinear optical processes in terms of interacting polaritons. Moreover 
equations describing the nonlinear interactions between polaritons become more similar to 
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those describing parametric interactions between photons widely adopted in quantum optics. 
Another, more fundamental reason, is that the standard second-order Born-Markov approx- 
imation scheme, usually adopted to describe the interaction with environment, is strongly 
bases-dependent, and using the eigenstates of the closed system provides more accurate re- 
sults. In order to obtain the dynamics for the polariton system we perform on the exciton 
and photon operators the unitary basis transformation 



= f/kSk; (10) 



being Pk = ( „ and 




= I I . (11) 

-^2k C2k 

In general photon operators obey Bose statistics, on the contrary the excitons do not posses 
a definite statistics (i.e. either bosonic or fermionic), but their behaviour may be well 
approximated by a bosonic-like statistics in the limit of low excitation densities. Indeed 

[Bn,Bl] = 5n',n-Y.'^l^^n'^ E ( I ^"^V^ • ^ + {N a \ dl^d_^ \ N (3)] I Na){Np \ . 

q N,a,l3 \ J 

(12) 

Thus, within a DCTS line of reasoning 23|], the expectation values of these transition oper- 
ators (i.e. I Na){Ni3 |) are at least of the second order in the incident light field, they are 
density-dependent contributions. Evidently all these consideration affect polariton statistics 
as well, being polariton linear combination of intracavity photons and excitons. As a conse- 
quence, even if polariton operators have no definite statistics, in the limit of low excitation 
intensites they obey approximately bosonic-like commutation rules. 
Diagonalizing Q,^: 

f/kOj^= = f2kf/k, (13) 

where 

\ u;2k 

uJi^2 are the eigenenergy (as a function of k) of the lower (1) and upper (2) polariton states. 
After simple algebra it is possible to obtain this relation for the Hopfield coefficients {2^: 

Xik = — C2k; C'lk = ^2k • (14) 



where 



X 



Ik 



c 



Ik 



V 



(15) 



V 



Introducing this transformation into Eq. ([9]), one obtains 



Pk = -^a^k + ^r-^^k''; 

where iV^^ = UTZ^^, which in exphcit form reads 



(16) 



Pik 

Ak 



-^^2kAk - ^S2k + ^2 k - ^R2^ ; 



(17a) 
(17b) 



where = ^cC^^^E'^, and = Xi^R^^, [m = 1,2). Such a diagonahzation is the 
necessary step when the eigenstates of the polariton system are to be used used as the 



starting states perturbed by the interaction with the environment degrees of freedom 
The nonhnear interaction written in terms of polariton operators reads 



3Q| 



^k ^ =E Al it) / 9r!i{t, t')P,^^ (t')P/k, {t')dt' : 



(18) 



being 



N, 



eff 



V 

nsat 



ci^^t - 1') + 



{v^5{t-t')-iF{t~t')^Xl^^ 
The shift Sk(t) is transformed into 

5™k(t) = Y.Pk,it) f_ {K^^5{t-t')-2zF{t-t') )P,^^{t')Pi^{t')dt' . 

and 



(19) 



ijl 



(20) 



1 



mk 



eff 



V 



+2\4x^ikp-^J"kp-^zk 



(21) 



Equation (I17p describes the coherent dynamics of a system of interacting cavity polari- 
tons. The nonhnear term drives the mixing between polariton modes with different in-plane 
wave vectors and possibly belonging to different branches. Of course there are nonlinear 
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optical processes involving modes of only one branch [8|,ll6|]. In this case it is possible to take 
into account only one of the two set of equations in (|T7|) and to eliminate the summation 
over the branch indexes in Eq. (1181) . 

Equations f[T7|) can be considered the starting point for the microscopic description 
of quantum optical effects in SMCs. They extend the usual semiclassical description of 
Coulomb interaction effects, in terms of a mean-field term plus a genuine non- instantaneous 
four-particle correlation, to quantum optical effects. Only the many-body electronic Hamil- 
tonian, the intracavity-photon Hamiltonian and the Hamiltonian describing their mutual 
interaction have been taken into account. The proper microscopic inclusion of losses through 
mirrors, decoherence and noise due to environment interactions will be the main subject of 
the following sections. 



III. QUANTUM LANGEVIN NOISE SOURCES : LAX THEOREM 

In order to model the quantum dynamics of the polariton system in the presence of losses 
and decoherence we exploit the microscopic quantum Heisenberg-Langevin approach. We 
choose it because of its easiness in manipulating operators differential equations, and above 
all, for its invaluable flexibility and strength in performing even multitime correlation calcu- 
lations, so important when dealing with quantum correlation properties of the emitted light. 
Moreover, as we shall see in the following, it enables, under certain assumptions, a (compu- 
tationally advantageous) decoupling of incoherent dynamics from parametric processes. 



In the standard well-known theory of quantum Langevin noise treatment [26|, |27|] greatly 
exploited in quantum optics, one uses a perturbative description and thanks to a Markov 
approximation gathers the damping as well as a term including the correlation of the sys- 
tem with the environment. The latter arises from the initial values of the bath operators, 
which are assumed to behave as noise sources of stochastic nature. Normally the model 
considered has the form of harmonic oscillators coupled linearly to a bosonic environment. 
The standard statistical viewpoint is easy understood: the unknown initial values of the 
bath operators are considered as responsible for fluctuations, and the most intuitive idea is 
to assume bosonic commutation relations for the Langevin noise sources because the bath 
is bosonic too. Most times these commutation relations are introduced phenomenologically 
with damping terms taken from experiments and/or from previous works. In other contexts 
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a microscopic calculation has been attempted using a quantum operator approach. Besides 
its valuable results as soon as one tries to set a microscopic calculation for interaction forms 

n 

different from a 2-body linear coupling [26f], e.g. acoustic-phonon interaction, some problems 
arise and one is lead to consider additional approximations in order to close the equations 
of motion and obtain damping and fluctuations. 

In 1966 Melvin Lax, with clear in mind the lesson of classical statistical mechanics of 
Brownian motion, extended the noise-source technique to quantum systems. In general, the 
model comprises a system of interest coupled to a reservoir (R). Considering a generic global 
(i.e. system+reservoirs) operator, a first partial trace over the reservoir degrees of freedom 
results in still a system operator, a subsequent trace over the systems degrees of freedom 
would give an expectation value. In order to be as clear as possible we shall denote the 
former operation on the environment by single brackets ( )r, whereas for the combination 
of the two (partial trace over the reservoir and subsequent partial trace over the system 
density matrices) the usual brackets ( ) is used. His philosophy was that the reservoir can be 
completely eliminated provided that frequency shift and dissipation induced by the reservoir 
interactions are incorporated into the mean equations of motion, and provided that suitable 
operator noise sources with the correct moments are added. In Ref. 28| he proposed for the 
first time that as soon as one is left with a closed set of equations of system operators for 
the mean motion (mean with respect to the reservoir) they can be promoted to equations 
for global bare operators (system+reservoir) provided to consider additive noise sources 
endowed by the proper statistics due to the system dynamics. He showed that in a Markovian 
environment these noise source operators must fulfill generalized Einstein equations which 
are a sort of time dependent non-equilibrium fluctuation-dissipation theorem. 

If a = {di, 02, ■ ■ ■} is a set of system operators, and 



{ M^) )r (22) 



dt 

are the correct equations for the mean, then one can show that the equations 

^ = i,(a) + F,(a,t) (23) 

are a valid set of equations of motion for the operators provided the additive noise operators 
F's to be endowed with the correct statistical properties to be determined for the motion 
itself. 
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The Langevin noise source operators are such that their expectation values {j^^)R vanish, 



but their second order moments do not 



28j. They are intimately linked up with the global 



dissipation and in a Markovian environment they take the form: 

( J-^(t) = 2{D^,)n 5{t - u) , (24) 

where the diffusion coefficients are 

Equation fl2^ is an (exact) time dependent Einstein equation representing a fluctuation- 
dissipation relation valid for nonequilibrium situations, it witnesses the fundamental cor- 
respondence between dissipation and noise in an open system. {D^u)r. becomes not only 
time-dependent, it is a system operator and can be seen as the extent to which the usual 
rules for differentiating a product is violated in a Markovian system. Equation and Eq. 
(!25|) make the resulting "fluctuation-dissipation" relations between D^^ and the reservoir 
contributions to be in precise agreement with those found by direct use of perturbation 
theory. This method, however, guarantees the commutation rules for the corresponding 
operators to be necessarily preserved in time. This result is more properly an exact, quanti- 
tative, theorem which gives relevant insights regarding the intertwined microscopic essence 
of damping and fluctuations in any open system. 

In order to be more specific, let us consider a single semiclassical pump feed resonantly 
exciting the lower polariton branch at a given wave vector kp. It is worth noticing, however, 
that the generalization to a many-classical-pumps settings is straightforward. The nonlinear 
term R^^ of Eq. (ITSl) couples pairs of wave vectors, let's say k, the signal, and k, = 2kp — k, 
the idler. A general result for quantum systems interacting with a Markovian environment 
is that after tracing over the bath degrees of freedom ones remains with system equations of 
motion in the bence of the environment plus additional phase-shifts (often neglected) and 



relaxation terms 
couple read 

d 



28j. The Heisenberg Eqs. f|T7|) . involving system operators, for the generic 



j^{PI)r = ^u^,{Pl)R + qIWrVI , (27) 
13 



where we changed shghtly the notation to underhne that pump polariton amphtudes Pkp are 
regarded as classical variables (C-numbers), while the generated signal and idler polaritons 
are regarded as true quantum variables. 

The nonlinear interaction terms in Eq. fl27|) reads 

" V 



—I 

9k 



C* +V^X* 

risat " 



XkXtX^^ . (28) 



Neff 

It accounts for a pump-induced blue-shift of the polariton resonances and a pump-induced 
parametric emission. In Eqs. ^I7\ only nonlinear terms arising from saturation and from the 
mean-field Coulomb interaction have been included. Correlation effects beyond mean-field 
introduce non-instantaneous nonlinear terms. They mainly determine an effective reduction 
of the mean-field interaction and an excitation induced dephasing. It has been shown 9| 
that both effects depends on the sum of the energies of the scattered polariton pairs. While 
the effective reduction can be taken into account simply modifying V^xj the proper inclusion 
of the excitation induced dephasing requires the explicit inclusion into the dynamics of four- 
particle states with their phonon-induced scattering and relaxation. In the following we will 
neglect this effect that is quite low at zero and even less at negative detuning on the lower 



polariton branch 29|]. The renormalized complex polariton dispersion tUk includes the effects 
of relaxation and pump-induced renormalization, cDk = t^k ~ iV^^"^^ /2 + \^kj\' , and 

^ / ^ * I |2 V ^ , |2 



Neff^risat " " risat 

+2Kx|Xk/|Xkn. (29) 

The damping term r^*"*"* here can be regarded as a result of a microscopic calculation in- 
cluding a thermal bath (see next Sect.). 

Following Lax's prescription we can promote Eqs. fl271) to global bare-operator equations 

|Pk = -^^kPk + QkPlVl^ + ^4 
jPl = iu^,Pl + glP^.Vl + J-pt^ . (30) 

However, in this form it is not a ready-to-use ingredient, indeed its implementation in 
calculating spectra and/or higher order correlators would be problematic because the noise 
commutation relations ask for the solution of the same (at best of an analogous) kinetic 
problem to be already at hand. This point can be very well explained as soon as one is 
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interested in calculating (Pj^Pk), i.e. the polariton occupation, where the mere calculation 
is self-explanatory. We shall need 

and the diffusion coefficient for the two operators in reverse order. Thanks to the structure 
above we can easily see that all the coherent contributions cancel out and only the incoherent 
ones are left. Anyway the important fact for the present purpose is that they are proportional 
to the polaritonic occupation, these coefficients will be explicitly calculated in Section |Vl 
The general solution of Eqs. (!30!) in the pump reference frame reads: 

P(t) = e/o^(*')^*'p(o) + f e^t'Mit")dt"j^^^,^ 

Jo 

m-\J'^'^ U-L^^^ I (32) 



M 



A*(k,r) oJ* 



p.: 



i2kp-k 




where 



l^2kp-k — — '^(l-^2kp~k — 2Ci;kp) 5 

the pump is Pkp = ^k ^ ^"^^^ 

p — pi p *^'^kp* 

2kp-k ~ 2kp-k'^ ' 

^ pt = > pt e ''p , 

^2kp-k ^2kp-k 

Ul' =Ui, + hy, |Pkp' 



vp I 



A(k,r)=^?kKp'- (33) 
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Equation (132|) can be written in a more explicit form by exploiting the following identity: 

e^'i' ^W^* = ai{tl, t2) [ ' A{t)dt + ao{tl,t2)I 

Jti 

A+(tl,t2)e^-(*i'*2) - A„(tl,t2)e^+("'*2) 



A+(tl,t2)-A_(tl,t2) 

ai{tl,t2) 



gA+(tl,t2) _ gA_(tl,t2) 

A+(tl,t2) - A_(tl,t2) 
A±(tl,t2) = /" A±(r)dr 

*2 / A_(tl,t2) 

diag[A(t)]rft = ^ ^ 
tl \ A+(tl,t2) 



= — , w = — . (34) 



Eq. ( ]32l) with Eq. (iMll provides an easy and general starting point for the calculation 

of multi-time correlation functions which are key-quantities in quantum optics. Taking the 
expectation values of the appropriate products it yields 

+ rfr|ci(r, t)p2(D^t pj(r)+^ rfr|c2(r, t)|^2(Dp^^t) , (35) 



here 



ci(ti,t2) = ai(tl,t2) / rfr( ^ -zJ^^) + ao{tl,t2) 

/•t2 

C2(ti,t2) = ai(tl,t2) / rfrA(k,r). (36) 



The two diffusion coefficients are proportional to the polariton occupation, i.e. we need 
as known input sources the very quantities we are about to calculate and a self-consistent 
solution seems unavoidable. Concluding, even if exact, Lax's theorem is of no immediate use 
for it simply rearranges the various ingredients to the microscopic dynamics in a different 
way. It seems worth noticing however that what up to now appears as a very formal and 
academic line of reasoning will be the clue for all the subsequent physical arguments ending 
up into an innovative approach to quantum optics in the strong coupling regime. Indeed, as 
we shall see in Sect. |Vl under certain assumptions we will be able to overcome the above 
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mentioned difficulty elaborating a (computationally advantageous) decoupling of incoherent 
dynamics from parametric processes. 

Anyway, it is the structure of Eq. fl31l) for the diffusion coefficients which allows, phys- 
ically speaking, to account for each contribution in its best proper way recognizing easily 
the dominant contribution. Indeed, reconsidering Eq. (PT]) in the light of the proper kinetic 
equation for the polariton population dynamics — the subject of the following section — , 
it is very clear that thanks to its structure all the coherent contributions cancel out auto- 
matically giving us an easy way to separate coherent and incoherent parts, but at the same 
time to treat them on an equal footing when calculating the final result. 

IV. MICROSCOPIC MARKOV CALCULATION OF POLARITON PHOTOLU- 
MINESCENCE 

Excitonic polaritons propagate in a complex interacting environment and contain real 
electronic excitations subject to scattering events and noise, mainly originating from the 
interaction with lattice vibrations, affecting quantum coherence and entanglement. For a 
realistic description of the physics in action, we need to build up a microscopic model taking 
into account on an equal footing nonlinear interactions, light quantization, cavity losses 
and polariton-phonon interaction. To be more specific as a dominant process for excitonic 
decoherence in resonant emission from QWs we shall consider acoustic-phonon scattering 
via deformation potential interaction, whereas we shall model the losses through the cavity 
mirrors within the quasi-mode approach (see Appendix J. It is worth pointing out that the 
approach we are proposing may be easily enriched by several other scattering mechanisms 
suitable for a refinement of the numerical results. 

In the view of the change of bases previously-mentioned, so imperative for a proper 
Markov calculation, we decide to treat the coupled system, described by the three Hamilto- 
nian terms He-, He and Hj, as our system of interest weakly interacting with the environment. 
In practice, this means to start from the linear part of the Heisenberg equations of motion 
in Eq. ([5]), which can be considered in the spirit of Sect. [TTllas system-operator equations, 
without the input term. Once obtained the polariton modes via a unitary diagonalizing 
transformation Eq. (JT3l) . we apply, to the coupling of this system with the environment, the 
usual many-body perturbative description. We end up with the customary Bogoliubov-Born- 
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Green-Kirkwood-Yvon (BBGKY) hierarchy which to the first order gives us the coherent 
input field, whereas to the second order the phonon and radiative scattering terms. As widely 
used in the literature, we shall limit ourselves up to this point thus performing a second-order 
Born-Markov description of the environment induced effects to the system dynamics. To 
exemplify our approach we shall calculate the relaxation rate of (O) = {P'l^Pa^ sole 
case of acoustic phonon interaction, any other scattering mechanism will be treated in the 
same way. The rate equation governing the incoherent dynamics to the lowest order of the 
polariton occupation is a relevant quantity we exploit in the next Sect, when we propose our 
DCTS-Langevin recipes for the calculation of multi-time many-body correlation functions. 
Being the full trace of the polariton density over the reservoir and the system degrees of 
freedom a relevant physical observable that we need to solve numerically, we prefer to give 
explicitly full account of the manipulations we have followed. It is worth underline that the 
very same formal treatment, i.e. Markov approximation, can be performed easily on the 



system operator arisen form the partial trace over the reservoir density matrix, ( )r 28|. 



In this latter guise damping and dephasing enter the mean system operator equation (l22ll . 



starting point for the Lax's theory of quantum noise [28||. A completely analogous procedure 
can be followed for the calculation of the dephasing rate of (Pak)- 

In the following the DOTS description of the interaction with the environment is limited 
to the lowest order. This means that effects like final state stimulation of scattering events 
are neglected. At the lowest order, the acoustic phonon interaction Hamiltonian can be 
expressed only in terms of excitonic operators as jsO] 

Kc-pH= E I 15 k + qll) (15 k I j (b^ + bl^ 

= J2Q^F^ + Q^F1^, (37) 
q 

is described in the Eq. (lA.Sp of the Appendix. Simbolically stands for all the system 

Hamiltonians, i.e. free dynamics and parametric scattering. The standard microscopic 



perturbative calculation 
-iO) = { - 



28| gives: 



X ( [O, Qq ] g-q(-n) ) - ( Q-^{-u) [O, ] ) , (38) 
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where the meaning of the new symbols are self-evident. 

Within the strong coupling region, the dressing carried by the nonperturbative coupling 
between excitons and cavity photons highly affects the scattering and for a microscopic 
calculation we are urged to leave the couple picture of Eqs. and move our steps into the 
polaritonic operator bases. Our aim is to produce a microscopic description of damping and 
fluctuation and to apply it in experiments with low and-or moderate excitation intensities, 
thus we expect the strong-coupling regime to become crucial in the scattering rates mainly 
through the polaritonic spectrum. In the spirit of the DCTS we shall consider them as 
transitions over the polaritonic bases obtained form excitons and cavity modes states and 
the linear diagonalizing transformation Eq. (IT^ . form exciton and photon operators to 
polariton operators, can then be rewritten as 

I 15k')(15k 1= 5^X,k'X;k I ^k')(jk I , (39) 
hi 

where it is understood we have transition operators on the left-hand side representing ex- 
citons, whereas on the right-hand side polaritons. Within the Born-Markov description we 
are left with 



dt 
with 



exc—ph tK 



El 4'' 



Qz 1 2 



1 , 1 



w, 



(sk),(rk') 



5(esk - erk' ± ^(k'-k,g,)) ( '^(k'-k.g,) + 2^2 



± 

(sk),(rk') 



± 

(rk'),(sk) • 



(41) 



rk' 



These happen to be the same ingredients used in Ref. [30| studying bottleneck effect in 
relaxation and photoluminescence of microcavity polaritons within a bosonic Boltzmann 
approach. 

Within the quasi-mode approach, the emitted light is proportional to the intracavity 
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photon number {tc is the transmission coefficient) : 

I^\t) = t^(4ak)(t) = tlJ2\C^^.\'{P!A){t) 



(42) 



By applying the whole machinery, when including only the lowest order terms in the input 
light field, the following equation for the polariton-occupations dynamics is obtained: 

d 



dt 



(43) 



Ik' 

with the generation rate given by 

9i,k= ^cCkp-£^k"*'^^'^k,kp(-Pi|k)co/i + ^cCkp-^k"^ ^'^k,kp(-Pj,k)co/i 

The phonon-emission (+) and phonon-absorbtion (— ) scattering rates read 



1 


k'-k 


f + fe°)2 


fmS \] 







(44) 



(45) 



X 



the 3D phonon wave vector is (q, q^) , whereas is calculated so that the energy conservation 
delta function 5{fux}jk' — hwik ± E^^) is satisfied, 

S = [DJ^{q,)l!{k' - k) - DX{q,)ll{k' - k)) , (46) 

with the overlap integrals 

-1 -3/2 



1 + 



2{me + ruh) 

dz\xeik){z)\'e"^^' 



qjox 2 



(47) 



We shall treat the cavity field in the quasi-mode approximation, that is to say we shall 
quantize the field as the mirror were perfect and subsequently we shall couple the cavity 
with a statistical reservoir of a continuum of external modes. This way on an equal footing 
we shall provide the input coherent driving mechanism (at first order in the interaction) and 
the radiative damping channel (within a second order Born-Markov description). 
The escape rate through the two mirrors {I =left, r =right) is 



(c) 

'^i,kll 



27r 



Q.kiiT Xj/ ^'^^KiTM ^l^skii('^)r> 

s=l,r 



(48) 
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and the corresponding noise term reads 

rSii =y|C.kiirE/ du;n'^{u;)6{Ql^[u;]-uj,,)h\gM^)\\ (49) 

s=l,r 

V. DCTS-NONEQUILIBRIUN QUANTUM LANGEVIN APPROACH TO PARA- 
METRIC EMISSION 

As already discussed, Eq. (15^ with Eq. provides an easy and general starting point 
for the calculation of multi-time correlation functions which are key-quantities in quantum 
optics. Thus it would be easy-tempting to wonder if, through some appropriate, thought 
and motivated physical considerations, we were given the noise sources as known inputs. 
Thanks to the structure of Eq. fl3Tl) for the diffusion coefficients we are allowed, physically 
speaking, to recognize properly the dominant contribution. In order to be more specific let 
us fix our attention on the explicit form of 2(£)ptp ): 

2(DptpJ(t) = Yl w^k,k'(PiPk')(t) + ■ (50) 

k' 

Inspecting Eq. (1431) . it results that in the low and intermediate excitation regime the main 
incoherent contribution to the dynamics is the PL the pump produces by itself, the effects 
on the PL of subsequent pump-induced repopulation arising from the nonlinear parametric 
part is negligible, that is to say the occupancies of the couple signal-idler are at least one 
order of magnitude smaller than the pump occupancy. This means that in Eq. fISU]) we can 
consider at the right hand side the solution in time of Eq. fl43p . i.e. only incoherent lowest 
order contributions. The other important diffusion coefficient reads: 

2{Dp^p^)it) = J2 w^k,k'(Pk'Pi)(t) + + 7^ = 

k' 

= ^k,k' ( {PlA')it) + 1 ) + + 7^ (51) 
k' 

We decide to use a sort of bosonic-like commutation relation in the equation above but only 
in the present situation, restricted only to this precise case and to the noisy background 
(PkPj^)(0), responsible for spontaneous parametric emission. The reason is many-fold. The 
two terms, Eq. flSlj) and the above noise background, will enter in Eq. fl35p multiplied by 
|c2p which contains the pump already twice. As a consequence their contribution must be 
to hnear order. Besides, Eq. (j43l) can be considered as the very low density limit of the 
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rate equation obtained from the picture of polaritons as bosons as obtained in Ref. |30 |. 
It witnesses that when the focus is devoted to the sole incoherent lowest order dynamics, 
bosonic commutation rules for polaritons may be employed, though carefully. Moreover, 
direct computation for normal incidence gives 

bI] = 5n',n - J2 '^nq'^n'q 5Z ( I ^i^q I + I ^-q^-'i I ^/^) ) I (^/^ I ' 

q N,a,f3 \ J 

(52) 

where as usual | NjS) and | Na) are N-pair eh pairs. Thus, within a DOTS analysis, the 
dominant term to the lowest order in the commutator is (5-like, whereas the two additional 
contributions, being proportional to the electron and hole densities, are nonlinear higher 
order corrections, contributing to the lowest order nonlinear dynamics but negligible for very 
low density, i.e. linear order. In the following we shall indicate as (-Pk-Pk)pL this solution 
representing the (incoherent) polariton occupation of the pump-induced PL. In the following 
subsection we show that this choice guarantees consistency between the rate equation Eq. 

and the complete solution we are about to present in Eq. (jSlD in the limit of pump 
intensity tending to zero, i.e. when it is the incoherent PL which governs the dynamics. 

polariton occupation dynamics 

With the notation introduced so far, the Heseinberg-Langevin equation governing the 
dynamics are those of Eq. fl30l) which we report here again: 



dt 

j/l = ^^kpI + 9lAvl + -^p; ' (53) 

where u^ = u^- iV^'"'^ /2 + K AX ' ^^^h 1^°*^ = (Pif'^) + 7^'^)- 
The general solution for the polariton occupation reads 
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) 



(54) 



In all the situations under investigation, the thermal population of photons at optical fre- 

(c) 

quencies are negligible, hence Fj^. ~ 0. Moreover, in the limit of pump intensity tending to 
zero it is the PL which governs the dynamics. Indeed Eq. (H3l) in this situation reads 



when, at least formally, we consider the right-hand side as known we can integrate, obtaining 



which is the limit of excitation intensity to zero of Eq. (IMI) . The form of 2 {Dp^p^ ) guarantees 
this fact for the reverse order calculation. 

In Eq. (1541) it is evident the great flexibility of the Langevin method, even in single- 
time correlations. It represents a clear way to "decouple" the incoherent and the coherent 
dynamics in an easy and controllable fashion. In the important case of steady-state, where 
the standard Langevin theory could at least in principle be applied, we have nonequilibrium 
Langevin sources which become: 



d_ 

dt 



k' 




= -rt'\PlP^) + E w^M'iPlA') + 



k' 



k' 



= -Tt'\P^Pl) + E W^k,k'(PkP;) + + 7k 



k' 



k' 



giving 



ri:°*)((ptp^)(t) + i). 
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i.e. the standard statistical viewpoint is recovered in steady-state. 

Concluding this section, standard Langevin theory gives some problems in dealing with 
interaction forms more complicated than the standard linear two-body coupling and some 
additional approximations are needed. Lax technique, on the contrary, provides us with the 
correct Langevin noise sources in the generic nonequilibrium case no matter of the operatorial 
form of the reservoir (weak) interaction Hamiltonian to implement. They properly recover 
the well-known steady-state result even if they depend, in the generic case, on the scattering 
rates rather than on dampings, on the contrary to standard Langevin description. 



time-integrated spectrum 

The spectrum of a general light field has always been of great interest in understanding 
the physical properties of light. Any spectral measurement is made by inserting a frequency- 
sensitive device, usually a tunable linear filter, in front of the detector. What is generally 
called "spectrum" of light is just an appropriatelx^jwrmalized record of the detected signal 
as a function of the frequency setting of the filter\fil\. Here we are interested to the power 
spectrum of a quantum-field originating from pulsed excitation and thus not at steady state. 



The time- integrated spectrum of light for a quantum-field can be expressed as |3l|, |32 |. 
Xk(^,r) = f dt, rrft2(4-^(ti)4^Ht2))e-(^-^-)(^-*^)e-(r+^-)(^-*^) , (55) 

J - J to J to 

where F is the bandwidth of the spectrometer (e.g. of the Fabry-Perot interferometer) and 
{E''^^) are the field operators corresponding to the light impinging on the detector, k 
is nothing but a proportional factor depending on the detector parameters and efficiency. 
Within the quasi-mode approach [ssl the spectrum of transmitted light is proportional to 
the spectrum of the intracavity field. In our situation, in the very narrow bandwidth limit 
and considering a beam with given in-plane wave vector, [34] it reads 

Jk(^,T) = r dt, r dt,{dl{t^)d^{t,))e-'^^'^-'^^ . (56) 

J- — h J to J to 

By expressing the cavity-photon operator in terms of polariton operators, one obtains 

Uu;,T) = r dt^ Trfta V |ak|'(i^Uii)Ak(t2))e-^"(*^-*^) . (57) 

J - J to J to 
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In our experimental conditions the upper polariton contribution is negligible, thus we need 
to calculate 

Jk(a;,T) = ^L^p9^ dt, f dh{Pl{t,)P^{h))e-^-^'^~'-^ . (58) 

^ - to J to J to 

By using Eq. (!32l) and the properties of noise operators (^^, one obtains: 

(P,I(tl)Pk(t2)) = Ci(0,ti)*Ci(0,t2)iVk(0) + C2(0,ti)*C2(0,t2)(iV2k,-k(0) + l) + 
rta 

/ dT ci(r,ti)*ci(r,t2) V Vrk,k'(^iA')pL(^) + (59) 
+St^Mnit^M) ^ dr C2(r, ti)*C2(r, W^2k,-k,k' ((PiPk')pL(r) + l) + 72k,- 

VI. NUMERICAL RESULTS 



In order to perform numerical calculations, we need to discretize in k-space. Although, 
thanks to confinement, cavity photons acquire a mass, it is about 4 order of magnitude 
smaller than the typical exciton mass, thus the polariton splitting results in a very steep 
energy dependence on the in-plane wave vector near k = (A; = ujsm6/c) (see Fig. 1). 
This very strong variation of the polariton effective mass with momentum makes difficult the 
numerical integration of the polariton PL rate-equations fl43p . For example if PL originates 
from a pump beam set at the magic- angle (see Fig. 1) or beyond, a small temperature of 5 K 
is sufficient to enable scattering processes towards states at quite higher fc- vectors, thus it is 
necessary to include a computational window in k-space, significantly beyond kp^mp- Usually, 
in finite volume numerical calculations, the k-space mesh is chosen uniform, but a dense grid 
suitable for the strong coupling region would result in a grid of prohibitively large number 
of points for (e.g. thermally activated) higher k, on the other hand a mesh well-suited for 
polaritons at higher k-values would consist of so few points close to k = to spoil the results 
gathered from the numerical code completely of their physical significance. Following Ref. 



35l | we choose a uniformly spaced grid in energy which results in an adaptive k-grid (in 
modulus), in addition, thanks to the rotating symmetry of the dispersion curve, we choose 
a uniformly distributed mesh in the angle 6 so that k = {k,6). Unfortunately, even if this 
choice allows for a numerical integration of the polariton PL rate equations (1431) . it provides 
an unbearable poor description of the parametric processes ( 1541) . The incoherent scattering 
events and the PL emission rates are strongly dependent on the energy of the involved 
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polariton states. On the contrary parametric emission, being resonant when total momentum 
is conserved, depends strongly on both the zenithal 9 and azimuthal angles which become 
poorly described by such adaptive mesh when the dispersion curve becomes less steep. 
Our DCTS-Langevin method enables the (computationally advantageous) decoupling of 
incoherent dynamics from parametric processes allowing us to make the proper choices for 
the two contributions whenever needed. 

In particular, we seed the system at a specific k and first of all we calculate the pump- 
induced PL by means of Eq. (H3|) . Because of the very steep dispersion curve and the large 
portion of k-space to be taken into account, in the numerical solution we need to exploit the 
adaptive grid above mentioned. Afterwards we use this pump-induced PL, (-Pk-^k)pL, as a 
known input source in Eq. flM|) where it is largely more useful to discretize uniformly in k. 



We consider a SMC analogous to that of Refs. [12|, [16 1 consisting of a 25 
nm GaAs/Alo.sGao.rAs single quantum well placed in the center of a A cavity with 
AlAs/Alo.isGao.ssAs Bragg reflectors. The lower polariton dispersion curve is shown in 
Fig. 1. The simulations are performed at T = 5K and the measured cavity linewidth is 
^7c = 0.26 meV. The laser pump is modeled as a single Gaussian-shaped impulse of FWHM 
r = 1 ps exciting a definite wave vector kp and centered at t = 4 ps. We pump with co- 
circularly polarized light exciting polaritons with the same polarization, the laser intensities 
/ are chosen as multiple of Jq corresponding to a photon flux of 21 /im~^ per pulse. We ob- 
serve that Ref. fl^ excites with a linearly polarized laser whose intensities are multiple of 
an Iq corresponding to a photon flux of 21 /im~^ per pulse too. In situations where the PSF 
and the MF terms dominate the nonlinear parametric interaction, there is no polarization 
mixing and two independent parametric processes take place, the first involving circularly 
polarized modes only and the second involving counter-circularly polarized modes. Thus 
for comparison with theory the effective density in those experiments is half the exciting 
density: / = /^/2. 

□ 

It has been theoretically shown[30|], that it is quite difficult to populate the polaritons in 
the strong coupling region by means of phono n-scattering due to a bottleneck effect, similar 
to that found in the bulk. Let us consider a pump beam resonantly exciting polaritons at 
about the magic-angle. Relaxation by one-phonon scattering events is effective when the 
energy difference of the involved polaritons do not exceeds 1 meV. When polariton states 
within this energy window get populated, they can relax by emitting a phonon to lower 
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energy levels or can emit radiatively. Owing to the reduced density of states of polari- 
tons and to the increasing of their photon-component at lower energy, radiative emission 
largely exceeds phonon scattering, hence inhibiting the occupation of the lowest polariton 



states. Actually t 
lar circumstances 



lis effect is experimentally observed only very partially and under particu- 



361 ] . This is mainly due to other more effective scattering mechanisms 371] 



usually present in SMCs. For example the presence of free electrons in the system deter- 
mines an efficient relaxation mechanism. Here we present results obtained including only 
phonon-scattering. Nevertheless the theoretical framework here developed can be extended 
to include quite naturally other enriching contributions that enhance non-radiative scatter- 
ing and specifically relaxation to polaritons at the lowest k-vectors[37|]. In order to avoid 
the resulting unrealistic low non-radiative scattering particularly evident at low excitation 
densities, we artificially double the acoustic-phonon scattering rates. However, acting this 
way, we obtain non-radiative relaxation rates that in the mean agree with experimental 
values. 

We now present the results of numerical solutions of Eq. (??) taking into account self- 
stimulation but neglecting the less relevant pump-induced renormalization of polariton en- 
ergies. 

Figures 2 and 3 show the calculated time dependent polariton mode-occupation of a 
signal-idler pair at k = (0,0) and at k = (2/^^,0), respectively, obtained for four different 
pump intensities in comparison with the time dependent pump-induced PL at the corre- 
sponding k. The pump beam is sent at the magic angle [4j {km — 1-44 ■ 10^ m~^) which is 
close to the inflection point of the energy dispersion curve and is resonant with the polariton 
magic angle is defined as the pump value needed for the eight-shaped curve 
of the resonant signal-idler pairs to intersect the minimum of the polariton dispersion curve. 
It is worth noting that the displayed results have no arbitrary units. We address realistic 
input excitations and we obtain quantitative outputs, indeed in Figs. 2 and 3 we show the 
calculated polariton occupation, i.e. the number of polaritons per mode. In our calculations 
no fitting parameter is needed, nor exploited (apart from the doubling of the phonon scat- 
tering rates). Moreover our results predict in good agreement with the experimental results 
of Ref. [l^ the pump intensity at which parametric scattering, superseding the pump-PL, 
becomes visible. It is clear the different pump-induced PL dynamics of the mode-occupation 
aX kx = 2 km (Fig. 3) with respect to that at the bottom of the dispersion curve of Fig. 
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2. Specifically a residual queue at high time values, due to the very low radiative decay 
of polaritons with k- vectors beyond the inflection point, can be observed. Furthermore we 
notice that already at moderate pump excitation intensities the parametric contribution 
dominates. It represents a clear evidence that we may device future practical experiments 
exploiting such a window where the detrimental pump-induced PL contribution is very low 
meanwhile we face a good amount of polaritons per mode. Indeed for photon-counting coin- 
cidence detections to become a good experimental mean of investigation we need a situation 
where accidental detector's clicks are fairly absent and where the probability of states with 
more than one photon is low. Our results clearly show that there is a practical experimental 
window where we would address a situation where all these conditions would be well fulfilled. 

We now focus our attention on the positive part of the ky = section at different pump 
powers. In Fig. 4 we observe the clear evidence of the build-up of the parametric emission 
taking over the pump-induced PL once the seed beam has become enough intense, in par- 
ticular we can set a threshold around / = 10 Jq (/^ = 20 Jq). As expected, the parametric 
process with the pump set at the magic angle enhances the specific signal-idler pair with 
the signal in = and the idler in kx = 2km- We can clearly see from the figure that at 
pump intensities higher than the threshold the idler peak becomes more and more visible 
for increasing power in agreement of what shown in Ref. [^] and Ref. 



. However, at so 



high kx values the photon component is very small and even if the polariton idler occupation 
is very high (as the inset if Fig. 4 shows), the outgoing idler light is so weak to give some 
difficulties in real experiments 12] . Moreover we can notice that the parametric process re- 
moves the phonon bottleneck in the region close to k = 0. An analogous situation occurs 
also in Ref. 36|, though with a different SMC, where it can be seen the bottleneck removal 



in k = due to the parametric emission. 

Fig. 5 shows the impact on the time integrated patterns of the calculated pump-induced 
PL. We consider, for different excitation intensities, the solutions of Eq. ( !54l) with and with- 
out the pump-induced PL occupations. As can be seen its inclusion does not result in an 
uniform noise background, but it seems to somewhat remember its incoherent nonuniform 
distribution (the one depicted in the corresponding curve, i.e. PLpump? in Fig- 4). As can be 
clearly gathered from the figure, the pump-induced PL has a non negligible contribution in 
a region in k-space resonant for the parametric processes. As a consequence at intermediate 
excitation intensities it adds up to the parametric part reaching a contribution even compa- 
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rable to the peak of emission set in k = O"*". Only beyond the above mentioned threshold the 
parametric emission is able to take over pump-induced PL and results in the great emission 
in the bottom of the dispersion curve of Ref. [41]. These results clearly show that PL emission 
does not become negligible at quite high excitation densities, but, being amplified by the 
parametric process, determines a redistribution of polariton emission displaying qualitative 
differences with respect to calculations neglecting PL. An interesting question regarding 
these phenomena could be related to the impact in the global spontaneous emission of the 
two contributions in Eq. namely that of the homogeneous part |c2(0,t)p and the one 

originating from noise operators in the time integral in the last line. In the inset of Fig. 5 
we have depicted the ratio of the homogenous solution with the global emission at k = 0^, 
calculated without {P^Pk)ph- Rather surprisingly when increasing the pump intensity the 
two contributions (homogeneous and particular) continue to have comparable weights, hence 
for a proper description of the spontaneous parametric emission they must be both included. 

The calculated time-integrated spectra of the outgoing light at k = (0, 0) obtained at six 
different pump intensities for an excitation at the magic angle km are shown in Fig. 4. It 
can be easily noticed a threshold around = 20 Jq in perfect agreement with the results 
in Fig. 4 and with Ref. [l^. For intensities lower than the threshold, the signal in = 
(with the corresponding idler in k^ = 2km) shows a quite large nearly Lorentian shape. As 
soon as the threshold is passed over, the spectrum starts to increase super-linearly with some 
spurious queues due to (calculated) asymmetric signal/idler damping values. Noticeably the 
spectra show an evident linewidth narrowing for increasing pump intensities witnessing the 
parametric emission build-up. For the sake of presentation in the inset we also present some 
normalized spectra which give immediate evidence of the build-up of a narrow linewidth 
beyond the mentioned threshold. 

VII. SUMMARY AND CONCLUSIONS 

Based on a DOTS theoretical framework for interacting polaritons (see Sect. 2), we have 
presented a general theoretical approach for the realistic investigation of polariton quantum 
correlations in the presence of coherent and incoherent interaction processes. The pro- 
posed theoretical framework combines the dynamics controlled truncation scheme with the 
nonequilibrium quantum Langevin approach to open systems. It provides an easy recipe to 
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calculate multi-time correlation functions which are key-quantities in quantum optics, but 
as shown here even for single-time quantities it provides a natural and advantageous de- 
coupling of incoherent dynamics from parametric processes. We have elaborated equations 
whose structure is analogous to those one obtains by means of bosonization {24]. However, 
thanks to the DCTS approach we have been able to obtain microscopically nonlinear co- 



efficients with great accuracy. In particular in Ref. [2J] the nonlinear coupling coefficient 
contains additional terms due to phase space filling providing an interaction strength larger 
of about a factor 3. We believe that the difference is due to the way how the Bosonization 
procedure has been adopted. As a first application of the proposed theoretical scheme, we 
have analyzed the build-up of polariton parametric emission in semiconductor microcavities 
including the infiuence of noise originating from phonon induced scattering. Our numerical 
results clearly show the importance of a proper microscopic analysis able to account for 
parametric emission and pump-induced PL on an equal footing in order to make quanti- 
tative comparison and propose future experiments, seeking and limiting all the unwanted 
detrimental contributions. Specifically, we have shown that already at moderate pump ex- 
citation intensities there are clear evidence that we may device future practical experiments 
exploiting existing situations where the detrimental pump-induced PL contribution is very 
low meanwhile we face a good amount of polaritons per mode. It represents an exciting and 
promising possibility for future coincidence experiments even in photon-counting regimes, 
vital for investigating nonclassical properties of the emitted light. 

APPENDIX: THE INTERACTIONS WITH RESERVOIRS 

A quasi-two-dimensional exciton state with total in-plane center of mass (CM) wave 



vector k may be represented as 39 1 



^''^)=7|E^"'''"'^^(^e - rlze,Zf,)al,^a,,,,\ 0) , (A.l) 



re.rh 



where f and 5* are the volume of the unit cell and the in-plane quantization surface, whereas 
^i/v r(^c/v,r) are creation (annihilation) operator of the conduction- or valence- band electron 
in the Wannier representation. Te/h = (rjl/^j, Ze/h) are to be considered coordinates of the 
direct lattice, |0) is the crystal ground state and R the exciton center of mass coordinate 
R = {nigYe + mhXh) I [p^e + ^ft.) with me and mh the effective electron and hole masses. The 
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Wannier envelope function Fx is normalized so that the integral over the whole quantization 
volume {V= S-L) of its square modulus is equal to 1. The lattice properties of GaAs-AlAs 
QW structures are in close proximity, thus the acoustic-phonons which interact with the 
quasi-two-dimensional exciton can be considered to have three-dimensional character. The 
electron-phonon interaction Hamiltonian resulting from the deformation potential coupling, 
can be written as 

HF-pH = E ( ^ ) (^ca;,k+q«c,k + ^.aU-fq«-.k) + bl^) . (A.2) 



k,q 



Here ajjy^ ^, ac/v,v. are creation and destruction operator of the conduction- valence- band 
electron in Bloch representation. Transforming from the Wannier to the Bloch representa- 
tion, we shall project (1A.2P into the excitonic bases. Moreover, since we are interested in 
the IS exciton sector A = (n, a) only 

Fis = Wis{rl - rl)xe{ze)Xh{zh) 



It yields 



W.sip) = ^2/{^)^M-p/<) ■ (A.3) 



k,k',g2 



4'"^= k|+gn _ _ DX{q,)ll{k' - k)) , (A.5) 



here 



2puV 

being /-*- and /" overlap integrals given in Eq. fH7j) . 

We treat the cavity field in the quasi-mode approximation, that is to say we shall quantize 
the field as the mirror were perfect and subsequently we shall couple the cavity with a 
reservoir of a continuum of external modes. The coupling of the electron system and the 
cavity modes is given in the usual rotating wave approximation 

Hgm= ih^jdujg^{uj)alE^~\uj,t) + H.c. , (A.6) 
k 

In passing form the air to the SMC, we change from a 3D to a 2D quantization, it means that 
in the coupling once either (k, kz) or (k, uo) is chosen the third follows consistently. We have 
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chosen the latter for simphcity in deahng with the Markov machinery. In the Hamiltonian 
g]i{uj) is the couphng coefficient, a sort of optical matrix element, E^^ \uj,t) and E^\uj,t) 
are the two propagating normal modes of the external light. Modeling the loss through 
the cavity mirrors within the quasi-mode picture means we are dealing with an ensemble of 
external modes, generally without a particular phase relation among themselves. An input 
light beam impinging on one of the two cavity mirrors is an external field as well and it 
must belong to the family of modes of the corresponding side (i.e. left or right). It will be 
nothing but the non zero expectation value of the (coherent) photon operator giving a non 
zero contribution on the 1*^* perturbative order. All the other incoherent bath modes will 
have their proper contribution in the 2"^^^ order calculations. 

It is worth noting that the treatment of the cavity losses as a scattering interaction is 
a result of the form chosen of the effective quasi-mode Hamiltonian. However, even if a 
model Hamiltonian, the quasi-mode description has given a lot of evidence as an accurate 
modeling tool and it is widely used in the literature. Let us call R the quasi-mode reservoir 
Hamiltonian. It can be shown that the ffist order (coherent) dynamics for a generic operator 
O under the influence of the coherent part of the quasi-mode ensemble reads, (see Eq. (1A.6I) ) 

='^T.Y.9pi 4'\^P't) )r ( [^^,4] ) + H.c. , (A.7) 

k p 

where g]i{uj) is the coupling coefficient, E^^ \i^,t) is the propagating normal mode of the 
external hght and Ylpdp^ \^p^t) )r is the superposition of all the possible coherent 
pump feeds. 

An interesting situation occurs within the assumption of a fiat quasi-mode spectrum, 
an approximation almost universally made in quantum optics |33|]. It makes Eq. (HHj) 
independent of the frequency: 

7S = E Yl^",kP^k.kP = |a,kP E 7l? , (A.8) 

i=l,r i=l,r 

where 7^'^'' is the (i-side) damping of the cavity without the quantum well. 
Thus 

5^7l? = 2vr5^kkP, (A.9) 

i=l,r i=l,r 

there are two situations: 
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equal damping: 'y^^^ = "fl"^ , we can define the transmission coefficient of the i-side 

1 2 ^i,k • ,2 



\9^M' = ^=<^ (A.10) 



we know the ratio: 



(m) ^ 

D '^rM I (m) 1 (m) 

^ J ^',k = TT7?7tot,k 

M ' (m) , (m) I - _^^(™) 

and the transmission coefficients follow. 

In the light of the definition of t^^, it becomes evident that the semiclassical coherent 
input feed could also be modeled from the beginning with an effective Hamiltonian: 

Hp = ^^5^(4"^ai - 4+)ak) , (A.ll) 

k 

where (the C— numbers) = ^^tc^p{ E^\Vtp,t) )^ represent the incoming coherent 

input beams j^. 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cam- 
bridge University Press, Cambridge, 2000). 

[2] V. M. Axt and T. Kuhn, Rep. Progr. Phys. 67 433-512 (2004). 

[3] B. Honerlage, A. Bivas, and Vu Duy Phach, Phys. Rev. Lett. 41, 49 (1978). 

[4] R. M. Stevenson et al., Phys. Rev. Lett. 85, 3680 (2000). 

[5] S. Savasta and R. Girlanda, Phys. Rev. Lett. 77 4736 (1996). 

[6] K. Victor, V. M. Axt, A. Stahl, Phys. Rev. B 51, 14164 (1995). 

[7] S. Savasta and R. Girlanda, Phys. Rev. B 59, 15409 (1999). 

[8] S. Savasta, G. Martino, and R. Girlanda, Solid State Communication, 111 495 (1999). 

[9] S. Savasta, O. Di Stefano, and R. Girlanda, Phys. Rev. Lett 90, 096403 (2003). 

[10] C. Weisbuch et al, Phys. Rev. Lett. 69, 3314 (1992); R. Houdre et al, Phys. Rev. Lett. 85, 
2793 (2000). 

[11] J. Erland et al, Phys. Rev. Lett. 86, 5791 (2001). 

[12] W. Langbein, Phys. Rev. B 70, 205301 (2004). 

33 



[13] P. Schwendimann, C. Ciuti, and A. Quattropani, Phys. Rev. B 68, 165324 (2003) 

[14] S. Kundermann, M. Saba, C. Ciuti, T. Guillet, U. Oesterle, J. L. Staehli, and B. Deveaud, 

Phys. Rev. Lett. 91, 107402 (2003). 
[15] C. Ciuti, Phys. Rev. B 69, 245304 (2004). 

[16] S. Savasta, O. D. Stefano, V. Savona, and W. Langbein, Phys. Rev. Lett. 94, 246401 (2005). 

[17] K. Edamatsu, G. Oohata, R. Shimizu, and T. Itoh, Nature 431, 167-170 (2004). 

[18] C. Diederichs, J. Tignon, G. Dasbach, C. Ciuti, A. Lematre, J. Bloch, Ph. Roussignol and C. 

Delalande, Nature 440, 904 (2006). 
[19] C. Ciuti, P. Schwendimann, and A. Quattropani, Phys. Rev. B 63, 041303 (2001). 
[20] B. Krummheuer, V. M. Axt, T. Kuhn, I. D'Amico, and F. Rossi, Phys. Rev. B 71, 235329 

(2005). 

[21] C. Ciuti, P. Schwendimann, B. Deveaud, A. Quattropani, Phys. Rev. B 62, R4825 (2000). 
[22] S. Portolan, S. Savasta, O. Di Stefano, F. Rossi, and R. Girlanda, larXiv:0706.224TV 2 [cond- 

mat.mes-hall], submitted. 
[23] Th. Osterich, K. Schonhammer, and L. J. Sham, Phys. Rev. Lett. 74, 4698 (1995), Phys. Rev. 

B 58, 12920 (1998). 

[24] C. Ciuti, P. Schwendimann, A. Quattropani, Semicond. Sci. Technol. 18, S279 (2003). 
[25] See e.g. D. F. Walls and G. J. Milburn Quantum Optics Springer- Verlag (Berlin-Heidelberg) 
1994. 

[26] L. Mandel, E. Wolf Optical Coherence and Quantum Optics (Cambridge University Press) 
1995. 

[27] G. W. Ford, J. T. Lewis, and R. F. OConnell, Phys. Rev. A 37, 4419 (1988). 
[28] M. Lax, Phys. Rev., 145 110 (1966). 

[29] S. Savasta, O. Di Stefano, and R. Girlanda, Phys. Rev. B 64, 073306 (2001). 
[30] F. Tassone,C. Piermarocchi, V. Savona, A. Quattropani, P. Schwendimann, Phys. Rev B. 56, 
7554 (1997). 

[31] J. H. Eberly and K. Wodkiewicz, J. Opt. Soc. Am. 67, 1252 (1977). 
[32] J. D. Cresser, Phys. Rep 94, 47 (1983). 

[33] C. W. Gardiner and M. J. Collett, Phys. Rev. A 31, 3761 (1985). 

[34] B. Renaud, R.M. Whitley and C.R. Stroud Jr., J. Phys. B 10, 19 (1977). 

[35] F. Tassona, Y. Yamamoto, Phys. Rev B 59, 10830 (1999). 

34 



[36] A. I. Tartakovskii et al, Phys. Rev. B 62, R2283 (2000). 

[37] G. Malpuech, A. Kavokin, A. Di Carlo, and J. J. Baumberg, Phys. Rev. B 65, 153310 (2002). 

[38] P.G. Savvidis et. al, Phys. Rev. B 62, R13278 (2000). 

[39] T. Takagahara, Phys. Rev. B, 31, 6552 (1985). 



35 



c 
o 

a. 
v> 

1- 

c 
d) 



1 .5220 
1.5215 
1 .521 
1 .5205 
1 .5200 
1 .51 95 
1.5190 
1 .51 85 
1.5180 







- 

^ si 








V / k 

1,1,1,1,1, 



1x10 2x10 3x10 4x10 

k (m-^) 



5x10 



6x10 



FIG. 1: (color online) Energy dispersion of the lowest polariton branch for the structure of Ref. [1/ 
consisting of a 25 nm GaAs/Alo.sGao.rAs single quantum well placed in the center of a A cavity with 
AlAs/Alo.isGao.ssAs Bragg reflectors. The pump at the magic angle and its parametric scattering 
(blue curve) are schematically depicted. The latter scatters two pump polaritons in a signal-idler 
couple at k = and k = 2km. The red curve symbolizes incoherent pump scattering at T = OK, 
e.g. due to acoustic phonon interaction. Because of the very steep dispersion curve in the strong 
coupling region due to strong coupling, for a pump beam set at the magic angle or beyond, even a 
small temperature of 5 K is sufficient to enable scattering processes towards states at quite higher 
/c-vectors, thus it becomes necessary to include a computational window in k-space, significantly 
beyond kp^^p making difficult numerical simulations. 
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FIG. 2: Calculated time dependent polariton mode-occupation at k = (0, 0) obtained at four 
different pump intensities in comparison with the time dependent pump-induced PL at the corre- 
sponding k. The laser pump is modeled as a single Gaussian-shaped impulse of FWHM r = 1 ps 
exciting a definite wave vector kp = {km, 0) and centered at i = 4ps. The calculated values are in 
good quantitative agreement with the measured value of Ref. [3l , with no fitting parameter needed 
nor exploited. 
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FIG. 3: Calculated time dependent polariton mode-occupation at k = (2km,0). It has been 
obtained under the same condition as Fig. 2. It is clear the different pump-induced PL dynamics 
with respect to Fig. [2l specifically a residual queue at high time values. 
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FIG. 4: Time-integrated outgoing photon emission intensity. The pump is set at kp = {km,0). It 
is clear the evidence of the build-up of the parametric emission taking over the pump-induced PL 
once the seed beam becomes higher than the threshold around = 20 Jq. Moreover the parametric 
process removes the phonon bottleneck in the region closed to k = 0. As expected the specific 
signal-idler parametric scattering with the signal in kx = and the idler in = 2km is favourite 
and at higher pump intensities dominates the light emission. The polariton idler occupations for 
some pump values are depicted in the inset. Although polariton occupation at kx = 2km is so high, 
its photonic component is very small resulting in a very weak outgoing light beam. 
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FIG. 5: The impact on the time integrated patterns of the calculated pump-induced PL for 
different excitation intensities is shown (the pump is set in kp = (kra,0))- Here the solutions of 
Eq. (j54p with and without the pump-induced PL occupations are depicted. On the contrary to 
what implicitly considered in previous phenomenological theories it results in a non-uniform noise 
background and hence its momentum distribution has to be included for a realistic microscopic 
calculation of the emission patterns. Moreover it is non- negligible in a region in k-space resonant 
for the parametric processes and hence at intermediate excitation intensities it adds up to the 
parametric part reaching a contribution even comparable to the peak of emission set in k = 0^ 
up to the threshold around = 20 /q • In the inset the ratio of the homogenous solution with the 
global emission at k = 0+ (both calculated without {Pl^Pk)pL) is depicted. For increasing pump 
intensities the two contributions (homogeneous and particular) in Eq. (|54|) still display comparable 
contributions, hence for a proper description of the spontaneous parametric emission they must be 
both included. 
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FIG. 6: Calculated time-integrated spectrum of the outgoing light at k = (0, 0) normalized with 
respect to the pump seed obtained at six different pump intensities for an excitation at the magic 
can be easily noticed a threshold around = 20 /q in perfect agreement with the 
results in Fig. 4 and with Ref. [l^. For intensities lower than the threshold, the signal in fc^ = 
(with the corresponding idler in = 2km) shows a quite large nearly Lorentian shape, as soon as 
the threshold is passed over the spectrum start to increase super-linearly meanwhile the linewidth 
decreases witnessing the parametric emission build-up. In the inset the normalized spectra at 
increasing pump powers (indicated by the arrow direction) are depicted evidencing even better 
the linewidth narrowing. We notice also some spurious queues due to (calculated) asymmetric 
signal/idler damping values. 
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